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Abstract. Nano and micro meter sized dust particles travelling through the 
heliosphere at several hundreds of km/s have been repeatedly detected by 
interplanetary spacecraft. When such fast moving dust particles hit a solid 
target in space, an expanding plasma cloud is formed through the vaporisation 
and ionisation of the dust particles itself and part of the target material at 
and near the impact point. Immediately after the impact the small and dense 
cloud is dominated by collisions and the expansion can be described by fluid 
equations. However, once the cloud has reached fim dimensions, the plasma 
may turn collisionless and a kinetic description is required to describe the 
subsequent expansion. In this paper we explore the late and possibly collisionless 
spherically symmetric unconstrained expansion of a single ionized ion-electron 
plasma using N-body simulations. Given the strong uncertainties concerning the 
early hydrodynamic expansion, we assume that at the time of the transition to the 
collisionless regime the cloud density and temperature are spatially uniform. We 
do also neglect the role of the ambient plasma. This is a reasonable assumption 
as long as the cloud density is substantially higher than the ambient plasma 
density. In the case of clouds generated by fast interplanetary dust grains hitting 
a solid target some lO'' electrons and ions are liberated and the in vacuum 
approximation is acceptable up to meter order cloud dimensions. As such a 
cloud can be estimated to become collisionless when its radius has reached 
fim order dimensions, both the collisionless approximation and the in vacuum 
approximation are expected to hold during a long lasting phase as the cloud 
grows by a factor 10®. With these assumptions, we find that the transition 
from the coUisional to the collisionless regime could occur when the electron 
Debye length within the cloud is much smaller than the cloud radius Ro, i.e. 
A = Ad/Rq <S 1- This implies a quasi- neutral expansion regime where the radial 
electron and ion density profiles are equal through most of the cloud except at 
the cloud-vacuum interface. The consequence of A being much smaller that unity 
implies that the electrostatic fields within a cloud generated by a dust impact 
on a neutral target is ~ 100 times weaker than in the case of grains hitting a 
spacecraft, where the positive potential of the target is strong enough to strip-off 
all the electrons from the expanding cloud leading to a "Coulomb explosion" like 
regime (e.g. Peano et al (2007) [l]). 
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1. Introduction 

The problem of the expansion of a plasma into vacuum has received much attention 
in recent years mainly in the context of the understanding of the expansion of plasma 
clouds generated by laser irradiated materials [U [31 [H [H [3]. The expansion of 
negatively charged dust particles in cometary tails [71 |H] and the expansion of the 
solar wind plasma into the wake region of inert objects such as asteroids or the 
moon [21 has also stimulated theoretical and numerical studies on the problem of the 
expansion of a plasma into vacuum. The impact of fast moving clusters of atoms or 
molecules on a solid surface are also known to produce expanding plasma clouds. 
In particular, dust particles, typically in the micro to nano meter range, hitting 
spacecraft at velocities up to hundreds of km/s have been repeatedly detected in space 
[TOl \TT\ fT2l \13\ [Ml [T5l [T6l fTTl fT8] . In most laser plasmas experiments only the electrons 
(not the ions) are heated by the laser's electromagnetic field. In these experiments, 
the initial state of the plasma is characterised by of hot electron population, carrying 
all of the energy, and a cold ion population, too tenuous for collisions to operate. 
On the other hand, dust impact generated plasma clouds are initially dense enough 
for electrons and ions to be in thermodynamic equilibrium everywhere within the 
cloud with the possible exception of the cloud/vacuum interface. One fluid models 
|191 1201 [m [52] , or two fluid models allowing for a separate description of ions and 
electrons (see the classical review paper [23]) are the appropriate tools to model the 
coUisional regime of the expansion. At some point however, provided the expansion 
takes place in vacuum or in a tenuous plasma and provided the coUisional mean free 
path of an electron grows faster than the cloud radius R, the expansion becomes 
coUisionless and a kinetic description necessary. Typically a nano meter dust particle 
impact at the solar wind speed is expected to ignite an expanding plasma cloud with 
some 10^ electrons and ions and a characteristic temperature of lOeV tmns coUisionless 
at cloud dimensions R > /.jm. The cloud then continues its expansion in the free 
coUisionless regime until its density has declined to a value comparable to the ambient 
plasma density, which at Earth's orbit occurs for a meter order cloud radius. Thus, 
during the free coUisionless expansion regime the cloud radius grows by a factor of 
order 10® before it merges with the ambient plasma. At such small scales magnetic 
forces can be safely neglected given that the Larmor radius of a low energy leV electron 
in interplanetary space already exceeds 600 m. 

The aim of this paper is to explore the unconstrained coUisionless expansion of a 
plasma cloud using kinetic simulations. The plasma made of identical single ionised 
ions and an equal number of electrons is assumed to be initially confined within 
a spherical shell which will be instantly removed to let the plasma expand freely 
into vacuum. The plasma is assumed to be initially at rest and at thermodynamic 
equilibrium implying equal ion and electron temperatures. This kind of initial 
condition differs from most of the published papers on the expansion of a plasma 
into vacuum where ions are generally assumed to be cold as in most laser heated 
laboratory plasma (e.g. [211 US] [201 [2Z] using fluid models, or [2H1 H] using kinetic 
models). We note that some of the presented results, in particular concerning the 
shape of the asymptotic electron velocity distribution function (see Figure [5]) , have 
been anticipated by Manfredi and Mola (1993) [29] using the strictly coUisionless 
Vlasov model and a hotter plasma. In this work we use a one dimensional N-body 
scheme to solve the equations of motion for a large number of ions and electrons so 
that collisions are not exclude a priori. The similarities between our results and the 
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results by Manfredi and Mola are then due to the fact that we assume the cloud to 
be marginally collisional from the beginning with the expansion further reducing the 
collisionality until an asymptotic "frozen" self-similar state is reached. It is worth 
noting that analytic self-similar solutions in the quasi- neutral limit, where charge 
separation is assumed to be small exist in the literature (see e.g. [30]). Such 
solutions are expected to hold in the limit of an electron Debye length being small with 
respect to the cloud dimension which we do suggest to be the case if the latter turns 
coUisionless during its expansion. However, the formation of an extended electron 
precursor, which appears to be inevitable in the spherical case |251 127j . suggest that 
the quasi-neutral assumption necessarily fails in the electron dominated, outer shells 
of the cloud. 

The plan of the paper is as follows. In Section[2]we introduce the basic parameters 
and equations relevant to the problem of the expansion of a spherically symmetric 
plasma. Since we focus our attention to the case of clouds which are initially collision 
dominated, we complete Section [5] with an outline the necessary conditions for a 
transition towards a coUisionless regime to occur. In section [3] the initial conditions 
and parameters for a simulation of a typical case are specified. The results of the 
simulation are discussed in detail in sections |4] to [71 A summary of the paper is 
presented in section [9] 

2. Definition of the problem 

We consider a spherically symmetric electron- ion plasma cloud made of N/2 singly 
charged ions and N/2 electrons expanding into vacuum. During the first phase 
of the expansion the plasma is supposed to be sufficiently dense to be dominated 
by interparticle collisions. This phase is conveniently described by fluid dynamics 
|19[ [211 [22l [23] and will not be treated in this paper. Under favourable conditions, 
however, the cloud radius R{t) may grow larger than the collisional mean free path 
of a thermal electron. At this particular radius R ^ Rq the cloud plasma becomes 
coUisionless and enters a new regime which is no longer tractable within the frame of 
a fluid theory. The reason for using electrons instead of ions to define the end of the 
collisional regime is that for a given temperature the collisionality of the latter may 
be substantially reduced as soon as the external shells of the cloud start moving faster 
than the ion thermal velocity as under such circumstances ions can no longer approach 
each other. As we shall see below the expansion velocity is indeed suprathermal for 
the ions after a short lapse of time roughly corresponding to the time required for the 
cloud to double its radius. 

2.1. Initial state of the cloud 

As the expansion is supposed to be collisional for R < Ro we assume that electrons 
and ions are initially (at time tg) in the state of thermodynamic equilibrium and 
confined within the spherical shell R = Rq. As the state of the plasma at the end 
of the fluid (collision dominated) phase is generally not known as it strongly depends 
on the cloud structure at the time of its formation and also on the adopted equation 
of state, we assume that ions and electrons are initially distributed according to a 
zero mean velocity Maxwell-Boltzniann distribution with temperature Tq. In order 
for the electron mean free path to be uniquely defined we do also assume that at t = 
the density within the cloud is spatially uniform. A more realistic description of the 
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initial state of the cloud when R~Rq, not even taking into account the presence of 
different ion species and neutrals (e.g. [21]), should include a non zero, radially varying 
fluid velocity profile, and a complex plasma- vacuum interface (see [121 )• We note that 
even if one assumes that the early, coUisional phase is governed by simple inviscid gas 
dynamics, the spatial and temporal structure of the expanding cloud has been shown 
to strongly depend on both the assumed energy equation (isentropic, isothermal, ...) 
and on the density and temperature distribution within the cloud at the time of its 
formation [211 [HI [22] . A few words on the radial velocity proflle u{r,t). In situations 
where the cloud radius R is allowed to grow much larger than the radius of the cloud at 
the time of its formation the evolution must be close to self-similar. Unlike the density 
and temperature profiles, the velocity profile then converges towards u = r/t for t ^ oo 
independently of the conditions at the time of formation and independently of whether 
the governing equations are fluid [U [1^1 [211122] or kinetic [Tl[151[17]- One may then be 
tempted to select a linear velocity profile u(r, t = 0) ~ r/St as initial condition for the 
collisionless regime, where —St is the instant of cloud formation. Unfortunately, such 
a profile is a priori incompatible with the assumption of a constant density profile 
unless very special, and unlikely, conditions exist at t = —St. More sophisticated 
initial conditions based on approximate self-similar solutions from compressible gas 
dynamics (e.g. [T9l[22]) will be discussed in a future publication. 



2.2. Basic parameters 

Previous works [25l [1] [27] on the spherical expansion of a plasma into vacuum 
have pointed out that the problem is characterized by the dimensionless parameter 
A = Xu/Ro at t = 0, where Ad is the electron Debye length (SI units) 

AD=f^)'''. (1) 

In equation ([!]) e is the absolute value of the electron charge, Eq the permittivity of 
free space, T the temperature and n the electron density. In equation ([!} and during 
the remaining of this paper we assume that temperatures are given in energy units, i.e. 
temperatures are implicitly multiplied by the Boltzmann constant fee- In situations 
when A <C 1 the thermal energy of the electrons is too low to allow for a substantial 
charge separation at the cloud surface: the expansion is quasi-neutral. On the other 
hand, in situations when A ^ 1 most of the electrons are energetic enough to overcome 
the electrostatic forces which bind them to the ions. In this case, the cloud becomes 
positively charged on a time scale of the order of Ro/vc {vc = (2T/mo)^/^ is the 
electron thermal velocity) and the associated peak electric field is much stronger than 
in the quasi-neutral case. In the limit A — > c« (the so called Coulomb explosion) all 
electrons escape from the cloud and the expansion is driven mainly by the repulsive 
forces pushing the unshielded ions away from each other. 

Let us now estimate the parameter A at time i = 0, when the plasma becomes 
collisionless. To this end we use the Fokker-Planck expression for the mean free path 
of a thermal electron 

= 16.4 (2) 

where A = ln(AD/?'s) is the Coulomb logarithm and Tg the strong interaction radius, 
usually defined as the larger of the classical distance for a strong electrostatic 
interaction between thermal electrons e^/(127reoT) or the de Broglie length for 
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a thermal electron fi,/(3moT)^/^, where h is the reduced Planck constant. For 
temperatures exceeding 9eV the quantum mechanical definition should therefore be 
used to define the rg. As typical cloud temperatures are expected to be of the order of 
a few eV up to at most 20eV and also because of the classical nature of the presented 
simulation, we stick to the classical definition throughout the paper. We emphasise 
that this assumption does not invalidate the subsequent discussions and the presented 
simulation for the case of temperatures higher than 9eV since we do only require rg 
being small with respect to the radius of the spherical shell rmin defining the inner 
boundary of the simulation domain (see . Now, even for an exceedingly hot cloud 
at 81eV, the quantum mechanical definition of rg is just three times larger than the 
classical definition. 

Equation Q is a good estimate of the mean free path of a thermal electron in 
a weakly coupled plasma with A > 10. For values 2 < A < 10 equation ^ may still 
be used as a fair estimate of Ic. By noting that the initial density uq is related to the 
cloud radius and the total number of electrons N/2 via 

»« - ^ (3) 

it follows from ([T|) and by setting lc{t = 0) = Rq in ^ that the dimensionless 
parameter A only depends on the total number of charged particles N in the cloud, 
viz. 

1/4 

(4) 



Rq 



A(7V) 



6iV 



We note indeed, that given the constraint Rq = lc{t = 0), the Coulomb logarithm 
A is a function of the total number of particles N via 6A^ = (4/3)"* e^'^A^'^ which 
leads to the relation A = 0.75Ae~^. Equation ([4]) indicates that the dimensionless 
parameter A is independent of the temperature Tq and much smaller than unity as N 
is generally a large number and A < 6 for < 10^. We conclude that at the time 
an initially coUisional plasma cloud becomes coUisionless it finds itself in the quasi- 
neutral expansion regime A ^ 1. For example, in the case of a 10~^" kg dust particle 
impacting on a spacecraft at solar wind speed the generated plasma cloud contains 
some N = 10^ charged particles corresponding to a Coulomb logarithm A ?» 5.5 and 
A = 0.017. 



2.3. From coUisional to coUisionless 

In this paper we restrict our discussion to the unconstrained expansion of a plasma 
cloud where the initial cloud's radius Rq increases by a large factor R{t)/Ro ^ 1 
before the dynamics of the expansion becomes affected by external factors, such as 
the ambient plasma. In the case of a dust impact generated plasma expanding into 
the ambient plasma (the interplanetary plasma) we can neglect the ambient plasma 
as long as the cloud density is much higher than the ambient plasma density. Indeed, 
for the nanodust impact considered above producing « 10^ charges with a typical 
per particle energy of lOeV (see [T7]) one finds, by setting 4 = -^o in equation 
that the plasma cloud can be considered coUisionless for Rq larger than a few ^m. 
Since the initial size of the cloud, just after impact, is necessarily comparable to the 
size of the dust grain (at most a few tens nm) most of the plasma within the cloud is 
coUisional at least during the first phase of the expansion. The growing cloud will then 
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turn collisioiiless provided the mean free path 1^ grows faster than the cloud radius R. 
Disregarding the weak dependence of the Coulomb logarithm A on the plasma density 
71 and assuming a polytropic equation of state T cx n''~^ one finds Ic/R (x rP^''~^/^\ 
Thus, for 7 < 4/3, the ratio lo/R is a growing function of R which means that the 
plasma must turn collisionless during expansion. If during the early phase of the 
expansion the conductive heat flux Q is unimportant the expansion must be adiabatic 
with an index 7a = 5/3 and the cloud turns increasingly collisional, at least as long 
as standard adiabatic fluid equations are applicable. Let us estimate under which 
conditions the conductive heat fiux is dominant by comparing the the conductive term 
and the adiabatic term in the energy equation for a spherically symmetric collisional 
gas of point particles: 

3DT T a , 2 ^ 1 a , 

ot;7 = — 2 7r^ " 2 7r^^ ^ 

2 Dt or n r^ or 

where D/Dt = d/dt + ud/dr is the material derivative. If the plasma in the cloud is 
collisional we can then use the Spitzer-Harm expression |33j for the conductive flux. 
Neglecting the contribution to the flux from the ions we therefore set 

Q = -1.6poWo (6) 

where pc = n^T is the electron pressure. Let us concentrate on the centre of the cloud 
at r = 0. If the plasma was initially uniform, at rest and spherically symmetric, it 
follows that density, pressure and temperature must have an extremum at r = and 
a vanishing first derivative at all times. On the other hand if the fluid velocity was 
initially zero at the centre it has to stay so for ever given that the acceleration is 
du/dr = —Q^^dpjdr = by virtue of the vanishing first derivative of the pressure. 
One may Taylor expand the velocity near r = as u{r) = M'(0)r + O(r^), where 
' = d/dr. According to the above discussion, the Taylor expansion of the temperature 
up to the first non constant term is T(r) = T{Q) + T"{0)r'^/2 + .... Of course, the 
same expansion can be applied to both density and pressure. To lowest order in r we 
can then write the energy equation ([5]) for the central region of the cloud as 

I ^ = -STou'iO) + 1.6vJ,T"{0) (7) 

where Tq = T{0). The first term on the right in equation Q corresponds to the 
adiabatic cooling due to the expansion while the second term corresponds to the non 
adiabatic cooling due to heat conduction. For the expansion to be non adiabatic 
the latter has to be of comparable order, or larger, than the former. In this case 
the effective polytropic index of the plasma is smaller than adiabatic 7 < 5/3 and a 
transition from collisional to non collisional becomes possible if conduction is strong 
enough to reduce 7 below the critical value 4/3. In order to estimate the relative 
importance of the 2 terms in (O we need an estimate of m'(0) and T"(0). Using the 
available macroscopic parameters of the cloud, like its radius R and the expansion 
velocity up of the front, it is natural to set u'(0) ~ uf/R and T"(0) « Tq/R^. We 
can then estimate the departure from adiabaticity by comparing the conductive to the 
adiabatic term, viz 

non adiabatic term 1.6 Wp L 
adiabatic term 3 up R 

From ([5]) it appears that the expansion is adiabatic in the limit of vanishing small 
mean free path le/R — > 0. However given that the typical expansion velocity up is of 
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the order of a few times the ion thermal velocity, the ratio Vc/u-p = 0([mi/mo]^/^) '3> 1 
is much larger than unity. Thus, unless the Knudsen number l^/R ^ 10~^ conduction 
is not negligible and the expansion cannot be adiabatic. Clouds with negligible heat 
conduction, i.e. for le/R <C 10~^ much smaller than (mo/mi)^''^, are expected both 
to remain collisional and to obey the equations of standard adiabatic gas dynamics 
in spherical geometry [34l [19] . This conclusion has to be softened somewhat as the 
collisional approximation for Q used in equation ([6]) may be inaccurate for Ic/R ^ 10"'^ 
(e.g. [35]) and is even expected to saturate at Qsat ~ 0.2pvc for Ic/R ^ 0.1 [361 137|. 



2.4- Problem reduction 

In the remaining of the paper we assume that the cloud under consideration goes 
through a phase where the Knudsen number Ic/R }^ 10^^ ensuring that the cloud 
turns coUisionless at some critical radius R = Rq. The free expansion then continues 
until the cloud's density has decreased to a level comparable to the ambient plasma. 
At Earth orbit, where the solar wind density is generally smaller than 10 electrons 
per cm"'^, the cloud density is substantially larger than the ambient plasma density 
for up to a meter order cloud radius R. Therefore, the expansion can be assumed to 
be both coUisionless and independent of the ambient plasma while the cloud radius 
R grows from Rq ^ fim up to i? ^ m, which corresponds to an expansion factor 
R/Rq = 0(10®). Given such a large expansion factor it is justified to assume that 
all particles within the cloud have purely radial velocities v = vf/r as the transverse 
velocity component v± (perpendicular to the radial direction) rapidly decreases during 
expansion since the angular momentum of individual particles L = mrv± is conserved 
in a coUisionless and spherically symmetric field. Neglecting the centrifugal force due 
to the transverse component of the particle velocity, the equations of motion for a 
particle of mass m and charge q in a spherically symmetric force field reduce to 

'i = ^Eir,t) (9) 
at m 

where E{r, t) is the radial electric field experienced by a particle at distance r from 
the cloud's centre. We shall verify a posteriori that neglecting the centrifugal term 
L^/ (m'^r^), which normally appears on the rhs of equation is justified by the fact 
that the field at the particle's position decreases asymptotically as r{t)^'^ (cf Section 
[6]), which is slower than the r{t)~'^ dependence from the centrifugal term. 

Given the spherically symmetric field E{r, t) assumed in ([9|) the particles must 
be interpreted as infinitely thin spherical shells rather than point particles. This 
approximation is justified as long as the number of particles within a given spherical 
shell is large with respect to unity, i.e. for radial distances ^ (3/4?™)^/'^ at time 
t = 0. Strict spherical symmetry reduces the original three dimensional system to a 
one dimensional system which can be treated much faster on a computer. The main 
drawback is an unrealistic description of the central part of the cloud which does not 
really matter as the small (and continuously decreasing) number of particles living 
in this region makes these particles statistically irrelevant anyway. Equations ^ and 
(jlOp must be supplemented by an equation for the electric field which for a distribution 
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E{r,t) = ^^, with Q(r,t)= ^ g,.. (11) 

ru(t)<r 

3. Setup and parameters of a selected simulation 




Figure 1. Snapshots of a subset of 3000 electrons in phase space. Initially, at time 
t = 0, the electrons are distributed uniformly within a spherical shell 0.1 < ^ < 1 
with radial velocities following a Maxwell-Boltzmann distribution with thermal 
velocity Vc = 8.34 (cf equation I I12I I. 



Figures [T] and [5] show the motion in phase space of a fraction of electrons and ions 
from a numerical simulation of a spherical cloud expanding into vacuum. Positions and 
velocities of N particles {N/2 electrons and N/2 ions) are time advanced according to 
equations ^ and (|10|) using a classical third order leap frog integration scheme |38| . 
The electric field is computed at every time step using the updated particles' positions 
and the field equation ([TT|) . 

The initial conditions consist in N = 80000 particles uniformly distributed within 
the spherical shell Tmin < r < Rq corresponding to A = 0.0538. Thus, even though 
the simulated number of particles N is much smaller than in a typical dust impact 
generated plasma cloud, it is still large enough for the key parameter A to be much 
smaller than unity so that the expansion is quasi-neutral. 

The inner sphere r < Tmin cannot be penetrated by particles and is merely there to 
avoid the divergence of the Coulomb potential for r when particles (actually thin 
spherical shells) approach the central region. In practice we chose rmin = 0.1i?0; which 
is both sufficiently small to minimise its influence on the overall system's evolution 
and sufficiently large with respect to the strong interaction radius rg to rule out binary 
collisions and self-charge effects. 
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Figure 2. Snapshots of a subset of 3000 ions in phase space. Initially, at time 
t = 0, the ions are distributed uniformly within a spherical shell 0.1 < ^ < 1 
with radial velocities following a Maxwell-Boltzmann distribution with thermal 
velocity vi = 0.834 (cf equation I ll2t . 



In the following, if not otherwise stated, we normalise charge to the elementary 
charge e, mass to the electron mass mp, length to r„iin, electric field to = 
e/(47r£:o'''min)' velocities to = e/(mo47reorjnin)^/^, time intervals to <„ = r^^^Jv-^ 
and temperatures to T,i = rricv'^. With these normalisation, and by consistently 
normalising density to r^^[^, the Debye length ([1]) reads Id = (r/47rri)^/^, the mean free 
path ^ Ic = T^/(7rnA) and the electric field of a point charge Q becomes E = Q/r^. 
We then set the initial temperature, for both electrons and ions, to Tq = 34.76 and 
the cloud radius to Rq ~ 10. The resulting Coulomb logarithm is then A = 4.03 and 
according to the mean free path 4 is equal to Rq = 10 as postulated. In code units 
the thermal velocity of the electrons is Ve = (2To)"'^/^ = 8.34 and the strong interaction 
radius ~ l/3To ~ 9.6 10""^ which, as required, is much smaller than both Rq = 10 
and r^in = 1. 

For convenience in Figure [T] and in all subsequent figures we use normalised 
positions ^ = r/R{t) with the temporal variation of the scale length defined by 
R{t) = i?o(l + i/io). We choose to set the arbitrary constant to = 10 in order 
to have dR/dt = 1. The ion to electron mass ratio is set to mi/rric = 100 so 
that ^0 does actually turns out to be of the order of the ion sound crossing time 
i?o/(3T/mi)^/^ = 9.8, a characteristic time for the initial system. As already stated, 
at t = particles are uniformly distributed within the spherical shell 0.1 < ^ < 1 
following Maxwell-Boltzmann velocity distributions for both ions and electrons: 

Mv) = ^e-(-M' (12) 
where Vj = (2To/mj)^/^ is the thermal velocity of the corresponding species j = {e, i}. 
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4. Asymptotic evolution, theoretical background 

The particle trajectories shown in Figures [1] and [2] illustrate two key aspects of the 
expansion which will be discussed in sections 14.11 and 14.21 First, as < — )• oo. all 
trajectories are seen to collapse towards the curve 11 = ^. As a consequence the 
temperature at a given position ^ is seen to decrease with time as the particles velocities 
appear to be less and less scattered as time progresses. Second, whereas ions rapidly 
line up in a structureless ribbon along the v ^ ^ curve, electrons converge towards a 
more complex structure, also aligned on the v — S, curve, but with a bulging of the 
ribbon in the region ^ < 2. As we shall see below the bulging is due to the bouncing 
motion of electrons trapped in an electrostatic potential well. 



4.1. Asymptotic convergence of particle trajectories 

In this section we show that for t — >■ 00 all particle trajectories must end up onto 
the V = curve provided the electric field decays sufficiently fast everywhere in the 
system. To this end we Taylor expand the asymptotic evolution of a particle velocity 
in terms of the small parameter ly = ti/t <C 1, i.e. v{iy) = v{0) + v{dv / dv)^^Q, where 
ti ^ to is just an arbitrary finite time level. From the equation of motion ^ we 
obtain the asymptotic evolution of a particle's velocity 



t^qE{^,t) 



t,\,^_mM (13) 

ti m ^ ,„ V r / m 

where Voc = w(t = 00). The asymptotic evolution of the particle's position is formally 
obtained by integrating (jlO|) . i.e. r(t) = 7*1 + J^^ v{T)dT. Using v{t) from equation 
(|13p. and for t ^ ti one obtains 

m-'i+v^-\fdr^-^^. (14) 
t t m 

The last term on the right hand side of equation (|14p vanishes for t ^ 00 provided E 
at position ^ decays faster than in which case = C(t 00) = Voo, confirming 
that the end point of a particle's trajectory lies on the — v curve. For a time 
dependence of the electric field £' cx it is possible to compute the slope of a particle's 
trajectory in the phase space directly from equations (|13p and (|14p . Indeed, assuming 
that for t 00 the variation of the electric field at a given particle position is due 
primarily to the time dependence of E rather than to the particle's motion, one obtains 
{v — Voa)/i£. — S.00) = 2 + a. Thus, for a = —2, corresponding to the final, self-similar, 
evolution of our system (see Figure[5]), (f — Woo)/(C~$oo) = 0, i.e. trajectories approach 
the S, ~ V curve on horizontal trajectories with v = const. In the particular case where 
E{£,,t) = Ei{S_)t\/t^ (which applies to the simulation for t/to ^ 1) equations (fTS]) and 
ifTt]) reduce to 

V = Voc — (15) 

in t 

,.j^ + !l_2M)|,„(l)^j_2M)|,„(l) ,10) 

t m t \ti J "I t \ti J 

Equation (fT6| shows that for sufficiently late times \£,{t)—(,oo\/£,oo ^ 1 confirming that 
the variation of the electric field at particle's position is asymptotically dominated by 
the field decay and not by the particle's motion. The interesting point about equation 
(|16p is that it shows that for t — >• 00 (which allows neglecting the ri/t term) particles 



Expansion of a spherical plasma cloud 



11 



approach their final position S^^o from the left or the right depending on the sign of qEi . 
Thus, in an overall positive electric field, which is indeed the case for the expanding 
cloud problem at hand (see Figure ^ , ions approach their final position from the left 
in {^,v) space while electrons approach their final position from the right. Thus, ions 
(electrons) which are initially on the right (left) of the v — ^ curve will first cross 
the V — ^ curve before converging towards their asymptotic position on horizontal 
V ~ const trajectories. This behaviour is already visible on the early phase of the 
expansion shown on figures [3] and |4l 



4-2. Shrinking of the volume occupied by particles in (^,u) space 

The shrinking of the phase space volume occupied by the particles in the (^, v) phase 
space is merely the consequence of the time dependence of the scaling length R{t). The 
equations of motion for an individual particle Q and (jlOp deriving from the general 
Hamiltonian of the system 

N 

H{ri,...,rN,Pi,-.,PN,t) = ^p]/2mj +qj(j){rj,t) (17) 

i=i 

where = rnjVj and —d(j)/drj ~ E{rj,t), it follows that any volume F = J dvdr must 
be conserved along particle trajectories in {v, r) space. Thus, F = R{t) J dvd£^ = const 
with the consequence that / dvd£, oc R{t)^^ = 1/t, i.e. the volume covered by the 
particles in the (^,v) phase space shrinks in time as 1/t. On the long term all 
particles must end up aligned on the v — ^ curve with the spatial distribution of the 
charges being a function of the initial conditions, i.e. on the dimensionless parameter 
A{N) = Xd/Ro only. 



5. Trapping, bouncing and freezing of particle trajectories 

Figures [3] and 2] show characteristic trajectories of selected electrons and ions in 
(^, v) phase space. From the figures it is immediately apparent that both species 
behave in a radically different way. Ions follow rather dull trajectories and are either 
accelerated outwards (in particular the outermost ones) or move at approximately 
constant velocity. The fastest electrons (in general the ones at largest radial distance ^) 
are seen to steadily reduce their outflow velocity in the attractive field of the positively 
charged interior of the cloud. However, electrons, with sufficiently low initial energy 
(the ones with end velocity v < 5), do cross the w = ^ curve and eventually bounce 
within an electrostatic trap. In order to understand the trajectories in the (^, v) space 
it may be useful to rewrite the equations of motion © and ([TU)) by setting v = v{£,{t)), 
viz 

e = ^ (19) 

where we have used dv/dt = £^dv/d£, with ^ = r/R{t) = r/t and R — 1. Equation 
shows that when a trajectory crosses the ^ = w it satisfies the condition dv/d£, ~ oo 
(cf Figures [3] and m unless the electric field E{£^,t) = 0. In this particular case the 
particle's velocity is constant v = Voo and shows that it takes an infinite time 
for the particle to reach the v = £, curve as ^(t) — ^(oo) oc 1/t. Multiple reflections 
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Electron trajectories 




^^0.0 0.5 1.0 1.5 

C = r/R{t) 

Figure 3. Sample trajectories of 102 electrons initially located in the range 
0.949 < 5 < 0.95 with velocities —12 < v < 12. The bottom panel is a zoom of 
the top panel. 

are associated with an equal number of crossings of the w = ^ curve, from top to 
bottom for an inward directed force and from bottom to top for an outward directed 
force. Note that particles approaching the centre ^ = make an artificial reflection 
there as their radial velocity must change from negative to positive. Given that in 
the simulations particles are not allowed to approach the centre at a distance less 
than ^min(i) = O.IRq/ R{t) reflection does effectively occur at ^min instead. Electrons 
bouncing back and forth in an expanding potential well can be efficiently cooled 
by first order Fermi deceleration. Now, despite the fact that both, bouncing and 
non bouncing electrons lose kinetic energy, the phase space volume occupied by the 
two populations evolves differently in time. Thus, whereas the velocity difference 
Av = I W2 — I'll between two electrons with initial velocities V2 and vi increases in time 
for vi,V2 ^ 10, the opposite is true for the trapped (and eventually bouncing) electrons 
with initial velocities vi,V2 ^ 10 (see Figure [3]). Loosely speaking, trapped electrons 
contribute to rise the particle concentration fj/Av in velocity space while non trapped 
electrons, and basically all ions, contribute to reduce the particle concentration in 
velocity space. Both effects are visible in Figure [5] which shows that fc{v) substantially 
grows over its initial value for v < 1.5 where Fermi deceleration occurs. For v > 1.5 
the electron trajectories in velocity space are divergent and fc{v) falls below its initial 
value, instead. 
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Ion trajectories 




-1.5 I 

0.0 0.5 1.0 1.5 2.0 2.5 

^ = r/R{t) 



Figure 4. Sample trajectories of 116 ions initially located in the range 0.949 < 
S, < 0.95 with velocities -2 < v < 2. 




Figure 5. Electron velocity distribution function foiv) at the end of the 
simulation at t = lOOto- Plotted as a reference, the initial distribution of the 
absolute value of the radial velocities /c ( | f | ) ■ 



A schematic representation of the electron dynamics in phase space and the 
corresponding evolution of their velocity distribution function fo{t) summarising the 
discussion of this section is shown in Figure [S] 
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5.1. Freezing of the bouncing motion 

Figure [3] shows that some electrons have time to bounce several times during the 
simulation before their motion becomes frozen onto the w = ^ curve. Others, for 
example the ones with initial velocity w w 10 and final velocity near tj w 3 are 
unable to perform a full bounce period during the time of the simulation. The 
reason is that the bouncing period th for a particle of mass m and charge q moving 
in the potential Q/r of a point charge Q is of the order th ^ 2TT{rar^ / qQY/^ 
(Kepler's third law) so that tb oc ^3/2^^3/2 ^ ^3/2 gj-Qws faster than the expansion 
time toxp = t/v = S^R/v. Thus, as the bouncing time to expansion time ratio 
ib/ioxp 2'K{v^mRo/qQy/^^^/^{R/Roy/^ grows as t^/^ the bouncing motion for a 
given particle will become frozen as soon as t\^{t, ^,Q)/t exceeds a value of order unity. 
Electrons for which tb/^cxp ~ 1 may just be a able to bounce once but may still lose 
a significant fraction of their initial kinetic energy if their total energy (kinetic + 
potential) is small and negative. A particle hitting a wall moving at velocity if is 
slowed down by Sv ~ —1w or bv ~ 2{w — v) depending on whether its initial velocity v 
is in the range v > 2w or 2w > v > w, respectively. Given that the expansion velocity 
of the electrostatic walls is rather slow compared to the electron thermal velocity 
Ve = 8.3 (the peak of the electrostatic profile in Figure [8] at ^ ^ 2 moves outwards at 
a velocity < 2), 2 or 3 reflections are required to slow down a thermal electron below 
the critical velocity w* « 2 which separates bouncing and non bouncing electrons 
(see Figures [5] and |6]) . The maximum initial velocity v* of a trapped electron can be 
estimated by equating the bouncing period t^ and the expansion time texp = 

27r [ ^ J ^ ' 

Equation (|20p may appear rather useless, Q{£,) being an unknown function of ^. 
However, we expect the right hand side of equation (|20p to be largest after a short 
time of the order Ao/we = 0.091 (the inverse electron plasma frequency) from the 
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beginning of the simulation. Thus, an a priori estimate is possible if we assume that 
Q corresponds to the number of electrons in the outer Debye shell of the initial cloud, 
i.e. Q ss 47rADi?o"-o = 1-5A^A. In the present case N = 80000 and A = 0.0538 so that 
Q « 6456. With these numbers and by setting R = Rq = 10 and ^ = 1 in equation 
(|20p. one obtains v* ^ A which appears to be a fair estimate of the upper limit for the 
initial velocity of trapped electrons (see Figure [3]). 

5.2. Ion distribution 

The asymptotic ion velocity distribution is shown in Figure [71 The figure clearly 
shows that the ion distribution closely follows the electron distribution for tj < 1.7 
(roughly twice the ion thermal velocity (2ro/"^i)^^^ = 0.83) while at higher velocities, 
up to the velocity of the fastest ions in the simulation v « 2.5, the ion density is 
in excess over the electron density. On the other hand, for v < 0.5 the asymptotic 
distribution fi{v, oo) falls below its initial value fi(v, 0) while the opposite occurs for 
V > 0.5. In principle, given that the mobile electrons tend to escape from the cloud, 
all ions should be accelerated outwards by the positive charge Q{£,) of the remnant 
(see bottom panel of Figure ^ and its associated, outwards directed electric field 
(see Figure [5]) . The outwards directed force should produce a displacement towards 
larger radial velocities of the original velocity distribution /i(|t'|,0) with fj{\v\,t) — 
below some minimum velocity. The displacement of /j towards higher radial velocity 
is visible in [7] for the fastest ions with end velocities v > 0.5. However, no region with 
/i = is visible at low velocities, though. The reason is that the slow ions are kept in 
the inner part of the cloud by the inwards falling, Fermi decelerated, electrons. The 
coupling between ions and cold electrons is so efficient there that the electric field 
asymptotically goes to zero for <^ < 1.5 (see Figure [H]). Figure [7] also shows that up 
to w « 1.7 both electron and ion velocity distributions are conveniently fitted by a 
single Maxwellian distribution with thermal speed Vm = (2To/Tn)^/^ = 1.17 where 
m = 0.5(?7ii + 771c) = 50.5 is the average particle mass. If we assume that the electron 
distribution is a Maxwellian sharply cut at an upper velocity v* , such that the missing 
electrons are those in the outer Debye shell of the initial sphere 1.5AN = 6456, we 
obtain the estimate Erf(w*/u^) = 1 — 1.5A (where Erf(a;) = 27r~-^/^ ds cxp(— s^) 
is the error function) which gives v* = 1.44. Obviously, 1.5AA^ is an overestimate 
of the number of electrons in the electron precursor. From Figure IH we take that 
this number is roughly 5 times smaller than 1.5AA'^ which allows for a more realistic 
estimate Eifiv* /ufn) = 1 - 1.5A/5, i.e. v* = 1.99. 

6. Charge distribution and electric field 

The spatial and temporal structure of the electric field is shown in Figure [8l As 
expected, after a dynamic initial phase, when particles are close to their asymptotic 
position in the (^,w) space, the electric field intensity decays as E{£^,t) oc t~^. The 
spatial structure E(S,,oo) is characterised by a central region ^ ^ 1 where, apart from 
fluctuations due to the small number of particles , the field intensity is essentially zero, 
corresponding to the region where the electron and ion density are equal (see Figure[9|). 
For C ^ 1-6 the field intensity rapidly rises towards the maximum £'max(l+i/iQ)^ w 2.5 
at ^ = 2.1, followed by a gentler negative slope over a much larger spatial scale. The 
large scale is associated with the scale of the electron precursor, which is obviously 
forged by the initial electron velocity distribution, and scales as ^faii « VcIq/Rq ~ 8.34. 
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Figure 7. Ion and electron velocity distributions at the end of the 

simulation at t = lOOto- Also plotted are the initial distributions /fc(|t'|,0) and 
the Maxwellian N{Trih^)^^^^ cxp[—{v/vrn)'^] where vjn = 
the thermal velocity based on the initial temperature To 
of mass m = (mi + me)/2 = 50.5. 
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Figure 8. Temporal evolution of the electric field profile. For t > 30 the spatial 
profile is essentially frozen while the amplitude declines oc f^^. 



The shorter scale of the rising part of the field profile is forged by the ions and scales 
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Vito/Ro = 0.83. 
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Figure 9. Radial linear density for ions and electrons (top panel). The bottom 
panel shows the linear charge density (thick black curve and right scale) in units 
of the positive elementary charge |e|. The thin blue curve corresponds to the net 
charge Q{^) contained within the the spherical shell of radius 



The charge distribution at the end of the simulation at t = lOOio is shown in the 
bottom panel of Figure [9] where, again, Q(^) represents the charge within the sphere 
of radius ^. The maximum Qmax ~ 1292 is reached for ^ « 2.7. Thus Qmax ~ O.SA^A, 
corresponding to roughly 1/5 of the positive charge contained in the outermost Debye 
shell of the cloud a.t t = to- On the other hand the charge at the position of the electric 
field maximum is Q{£^ = 2.1) « 1053 « 0.2457VA giving a maximum field intensity 
-Einax(V*o)^ ~ 0.245iVA/(2.1i?o)^ ~ 2.4 as confirmed by the latest profiles of E shown 
in Figure [H Changing to SI units we obtain a more useful expression 

£;max,SI = 8 10-"^ (21) 

where -Emax.si is expressed in V/m and R in metres. 



7. Ion and electron fluid velocities 



Figure [TO] shows the spatial fluid velocity profiles Uc.p for both electron and ions at the 
end of the simulation. Not surprisingly both populations have velocities profiles close 
to the V = curve (top panel). Plotting the fluid velocities Uc^p — C versus ^ (bottom 
panel) shows that while all ions and electrons located within ^ < 1.7 do closely follow 
the V = curve (meaning that they are already frozen) electrons with 1.7 < ^ < 6 stay 
below the v — ^ curve. These electrons are still flowing (falling) inwards along the 
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Figure 10. Fluid velocities Uj{^) = v fj{v,i;)dv/nj^ for ions and electrons 
at the end of the simulation (top panel) and uj — ^ for at two different times 
(bottom panel). The grid of dashed lines in the bottom panel gives the trajectories 
V = const individual particles are expected to follow in case of an electric field 
declining as t~'^ (see Section BTTI l . 



dashed lines, corresponding to the v = const trajectories predicted for a t^^ decHning 
electric field (see Section |4TT|) . The electron inflow velocity at i = 100<o peaks near 
^ « 2.3 at about 10% of the absolute ion fluid velocity. The ion velocity will not 
evolve signiflcantly after t = lOOio as it is already well approximated by Ui = ^. The 
asymptotic position of the electrons near ^ Ri 2 at t = lOOto, will lie some 10% closer 
to the centre of the cloud. This late displacement of the electrons will not modify 
the final structure of the electric field signiflcantly as the density of these inflowing 
electron is substantially lower than the ion density near ^ = 2 (see Figure top panel 
in [9]). On the other hand, electrons at ^(lOOto) ^ 3 are too fast for the ions to catch 
up with and constitute the flnal electron precursor. 

8. Application to the case of clouds formed by interplanetary 
nanoparticles impacts 

When dust particles travelling in the interplanetary space hit a solid object at 
characteristic velocities of the order of tens to hundreds km/s, a plasma cloud is 
generated at the impact point. The cloud is formed due to the vaporisation and partial 
ionisation of the dust particle itself and the target material. The subsequent expansion 
of the cloud is hemispherical rather then spherical as assumed in the above model. In 
the case of a non conducting and charge neutral target the above results should not be 
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modified in any substantial way. On the other hand, one should be extremely cautious 
when trying to interpret the electric signals measured on a spacecraft as spacecraft are 
generally positively charged due to electrons being stripped from their metallic surface 
through photoionisation by solar radiation. The associated electric field (typically of 
the order of a few V/m) exceeds the cloud's internal electric field very early during 
the expansion, causing stripping of most of the cloud's electrons [17]. Not only is the 
expanding cloud subject to charging so that the long term expansion is more like a 
Coulomb explosion[TJ[27] rather than a quasi- neutral expansion as described above, but 
the role of the photoelectrons, continually emitted and recollected by the spacecraft, 
must be taken into account when trying to interpret the voltage pulses measured on 
spacecraft antennae. Indeed, in the scenario described in appendix A of the Zaslawsky 
et al paper |39j). the voltage pulses measured on individual antennae in conjunction 
with the impact of nano dusts on the STEREO spacecraft are not a direct measure of 
the field within the post impact expanding plasma cloud but are the consequence of the 
equilibrium photoelectron return current towards the part of the antenna within the 
plasma cloud being interrupted. The interruption of the photoelectron return current 
induces an accumulation of positive charges on the antenna which is then measured 
by on board detectors as a temporal variation of the potential between the antenna 
and the spacecraft. 
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Figure 11. Electrostatic potential V{i;) and total charge Q(§) as in figure[9] 

Let US estimate the electrostatic potential through a real cloud generated by a 
nanoparticle hitting a target in interplanetary space using the above model. As we 
shall see the intrinsic field of the cloud is too weak to account for the impact associated 
potential pulses observed on STEREO. 

The electrostatic potential through the simulated plasma cloud is shown in Figure 
[TTl The total drop in the electrostatic potential is Ay w 550 while, as already stated, 
the peak charge is Qmax ~ 1292 (see bottom panel in Figure [9]). Assuming a linear 
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relation between AV and Qmax we find 

AV{N) w « 0.128 NA (simulation units). (22) 

2.35 

This relation can be interpreted as the electrostatic potential at the surface of a sphere 
of radius 2.35i? enclosing a charge Qmax- Multiplication of Equation (p2)) by e/ (47reo^) 
in SI units, leads to the dimensional version of (22): 

AVsi(A^, R) w 1.84 10"^" ^ (23) 

Jx 

where R is in meters and Vgi in Volts. 

Let us use equation ([^^ to estimate the voltage pulse AV^ due to the impact of 
a md = 10"^" kg grain (size ~ 10 nm) travelling at w = 300km/s. The number of 
electrons and ions within the plasma cloud can be estimated via the semi-empirical 
formula N/2 « 0.7m\-°'^v^ '^ye (see [III EI]), where [e] = C, [md] =kg and [v] =km/s. 
With these parameters we obtain « 1.4 10^ and, from the relation between A and 
A'' established in the discussion following equation ^ we find A(1.4 10^) = 0.016. 
Plugging this value into in equation ^ leads to Aysi(l-4 10^) w 4.1 10"^/i?. Thus, 
when the cloud's radius has grown to i? = 50cm, i.e. when its density has fallen 
to a value comparable to the interplanetary density, AVsi ~ 0.082mV only. Noting 
that the effectively measured voltage is obtained by averaging over the whole antenna 
length, i.e. by multiplying the above voltage by l/L, where L is the length of the 
antenna (6 m on STEREO) and / the length of the part of the antenna within the 
cloud, the voltage predicted by the model is far too week to be directly detectable and 
in any case much too weak to induce the up to lOOmV pulses observed on STEREO. 
Previous estimates of the voltage pulse associated with nano-size dust impacts were 
based on the assumption that the charge separation within the expanding cloud is 
total with ~ N/2 electrons in the precursorHJ. The assumption, adopted here, 
that charge separation only occurs at the time when the expanding cloud becomes 
collisionless, without influence of potential external fields, reduces the number of 
electrons in the precursor to a much smaller number of order 0.3A'^A. Using equation 
pip we can estimate the maximum electric field within the above cloud to be of order 
Emax,si ~ 1-8 10~^i?~^. Thus, when the cloud has reached R = 4.2mm the maximum 
field is already down to IV/m, i.e. comparable to the spacecraft's own field. Beyond 
this size, the cloud will start loosing its electrons to the spacecraft. However, even 
in case of a complete charge separation such that the electrostatic field outside the 
cloud is given by the Coulomb field Q/r, the latter is much too weak (at most a few 
mV) to account for the on board measured voltage pulses associated with nanodust 
impacts. As we shall explain in more details in a forthcoming paper, but as already 
briefly exposed in the appendix A of |39j . the strong potential pulses measured on 
individual antennas on STEREO are due to the perturbative effect of the expanding 
cloud on the photoelectrons surrounding the antenna. 

9. Summary and conclusions 

We have explored numerically the unconstrained spherically symmetric expansion 
of an initially uniform, overall neutral and at thermodynamic equilibrium cloud of 
immobile plasma. The initial temperature and density of the plasma are such that 
the cloud's radius equals the Fokker-Planck coUisional mean free path of a thermal 
electron, representing an admittedly crude model of an expanding cloud at the time 
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it becomes collisionless. Consistently assuming that the ion and electron velocity 
distributions are Maxwellians at the time of the collisional to collisionless transition, 
it follows that the key parameter of the problem A = Ad / Rq can be written as a 
function of the total number of ions and electrons in the cloud only. Due to the 
A oc A^~^/^ dependence (see equation (jH)) typical nano size grain impacts which are 
expected to ignite plasma clouds with N = 0(10*), A is always much smaller then 
unity, i.e. the collisionless expansion is quasi-neutral. 

During the initial phase of the collisionless regime most electrons (the less 
energetic ones) lose nearly all of their kinetic energy thorough Fermi deceleration in 
the expanding potential (see Figure [TT|) . On the other hand the outermost ions near 
the electric field maximum (see Figure |8]), are accelerated outwards by the positive 
electric field. The net effect is that ions and electrons do asymptotically tend towards 
having the same velocity distribution up to a threshold of the order of the ion thermal 
velocity (see Figure [7]). The fact that all particle trajectories converge towards the 
V = ^ curve as t — ?> cxd implies that electron and ion fluid velocities end up being simple 
linear functions of the distance r from the expansion centre. 

At late times the ion density profile is conveniently described by ni{r) oc 
exp[— (r/twm)^] where vjn = (2To/m)^/^ is the thermal velocity based on the initial 
temperature To and the mean mass m = O.birm + me). Because of A <C 1 the electron 
density Ue closely follows the ion density up to a distance r* solution of the equation 
Erf (r*/t?>Yn) « 1 — 0.3AiV. Beyond this point the electron density is flat up to a 
radial distance of the order of v^t as observed in Vlasov simulations [29]. The electric 
field is essentially zero for r < r* (see Figure E]) but rises towards a maximum on 
the ion length scale rrise ^ Vit and decreases slowly on the electron precursor length 
scale rfaii ~ Vct. At late times the maximum field intensity -Bmax is essentially nailed 
down by the number of electrons in the outer shell of the cloud of thickness Ad- 
This number can be expressed in terms of the total number of particles N and the 
diniensionless parameter A = Xd/Rq to give l.SA^A, implying that the electric field 
intensity must be proportional to NA/R^. In our representative simulation, we find 
-Bmax.si ~ 8 10~^^NA/R^ which we expect to hold as long as TV ^ 1 and provided the 
cloud has changed from collisional to collisionless during expansion. The electrostatic 
potential through the cloud has been found to be AVsi{N, R) « 1.84 10"i°7VA/i?. 

The electrostatic potential difference between the cloud's centre and infinity 
predicted by the model is far too weak to be account for the voltage pulses, sometimes 
exceeding lOOmV, observed on the S/WAVES TDS detector on the STEREO 
spacecraft following a nano-sized dust particles impact. One plausible reason for this 
discrepancy is that spacecraft are positively charged due to photoelectron emission 
through their sunlight exposed surfaces. The resulting electric field, typically of the 
order of a few Vm~^ at lAU from the Sun, generally exceeds the maximum field 
intensity within the plasma cloud before its dilution in the ambient plasma. The late 
evolution of the cloud is therefore dominated by the spacecraft field which strips most 
or all of the electrons from the cloud which then sees both its charge and its internal 
electrostatic potential field increase by a factor of order A^^ ^ 1. However, even in 
case of an unrealistically large nanodust impact generated cloud with some N = 10* 
and with all electrons stripped-off, the total electrostatic potential difference would 
merely be a smallish 20mV at the time of its maximum extension, when R w Im. 
Considering that the measured field is down by at least a factor R/L, where L is the 
the total length of the antenna {L = 6m on the Stereo spacecraft), one must conclude 
that the on board measured fields are not a direct measure of the cloud's intrinsic 
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field. Indeed, recent findings by Zaslavsky et al [31] indicate that nanodust impact 
associated clouds strongly affect the photoelectron environment of the antenna. In 
the scenario proposed by Zaslavsky et al the photoelectrons emitted by the sunlight 
exposed surface of the antenna are temporarily hindered to fall back onto it because of 
the presence of the cloud's perturbing field. The resulting net photoelectron current is 
strong enough to allow for a fast positive charging of the antenna which is compatible 
with the measured field intensities. The bottom line is that the presented model is 
not directly applicable to the case of plasma clouds generated by nanodust impacts on 
spacecraft as neither the spacecraft potential nor the surrounding photoelectrons have 
been considered. The model is however expected to be applicable in case of nanodust 
impacts on uncharged targets. 
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